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ABSTRACT 

We use a large sample of 174 Wolf-Rayet (WR) galaxies drawn from the Sloan Digital Sky Survey to study 
whether and how the slope of the stellar initial mass function depends on metallicity. We calculate for each 
object its oxygen abundance according to which we divide our sample into four metallicity subsamples. For 
each subsample, we then measure three quantities: the equivalent width of H/3 emission line, the equivalent 
width of WR bump around 4650A, and the WR bump-to-H/3 intensity ratio, and compare to the predictions of 
the same quantities by evolutionary synthesis models of Schaerer & Vacca. Such comparisons lead to a clear 
dependence of the slope of initial mass function (a) on metallicity in that galaxies at higher metallicities tend 
to have steeper initial mass functions, with the slope index ranging from a ~1.00 for the lowest metallicity of 
Z = 0.001 to a ^3.30 for the highest metallicity Z = 0.02. We have carefully examined the possible sources 
of systematic error either in models or in our observational measurements and shown that these sources do not 
change this result. 

Subject headings: galaxies: abundances - galaxies: starburst - stars: Wolf-Rayet 



1. INTRODUCTION 

The stellar initial mass function (IMF) describes the relative 
frequency with which stars of various masses are formed, and 
so controls nearly all aspects of the evolution of stellar sys- 
tems. An understanding of the IMF has thus long been one of 
the major goals in many astronomical fields, from theories of 
star formation to the interpretation of integrated properties of 
galaxies at the highest redshift. 

The studies of IMF began with the publication of the first 
and still most famous paper on the subject by Salpeter (1955), 
who found that the solar neighborhood IMF can be approx- 
imated by a declining power law, dN/dM <x M~ a , with a 
slope of a = 2.35. From then on, there have been consid- 
erable observational studies aiming to determine the form of 
the IMF (e.g. Scalo 1986, 1998; Kroupa 2001, 2002; Reid 
et al. 2002; Chabrier 2003). These studies have established 
that, for masses above 1 M®, the IMF can generally be quan- 
tified using the Salpeter form, while for low masses it deviates 
from this form, by flattening below 0.5 M , peaking at a char- 
acteristic mass of ~ O.IM,.) and then significantly declining. 
However, despite much progress, there are still many prob- 
lems that are far from being solved (Kennicutt 1998; Bate & 
Bonnell 2005), out of which the universality or the variation 
of the IMF is among the most fundamental and is still a con- 
troversial matter. In particular, it is not very clear whether 
or not the IMF, including both its characteristic mass around 
O.1M and its power-law slope at higher masses, depends on 
environmental conditions such as metallicity. 

One reliable way to derive an IMF is to count individual 
stars when both photometry and spectroscopy are available. 
Massey et al. have applied this method to numerous young 
star clusters and OB associations in the Galaxy and the Mag- 
ellanic Clouds (Massey et al. 1995, Massey 1999). This ap- 
proach is, however, infeasible for deriving the starburst IMF 
because the distance of the closest starburst galaxies are ^ 
2Mpc, where the spectrum of individual stars cannot be ob- 
tained. Thus the stellar content of galaxies at large distances 
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can only be derived from their integrated spectra (Leitherer 
1998). In such studies, many authors have preferred to Wolf- 
Rayet (WR) galaxies, where broad emission lines from WR 
stars are observed in the integrated spectra. WR galaxies 
have long been known, beginning with the discovery of such 
features in the spectrum of the galaxy He 2-10 by Allen, 
Wright, & Goss (1976). Osterbrock & Cohen (1982), and 
Conti (1991) introduced the concept of WR galaxies. These 
galaxies are thought to be undergoing present or very recent 
star formation that produces massive stars evolving to the WR 
stage, and so are ideal objects for studying the early phases of 
starbursts, determining burst properties and constraining pa- 
rameters of the IMF (Schaerer et al. 1999; Guseva et al. 2000). 

Many studies have been done trying to constrain the main 
burst parameters, such as the age and duration of the bursts 
and the IMF, by comparing the observed WR features with 
that predicted by population synthesis models for young star- 
bursts (Contini et al. 1999; Schaerer & Vacca 1998, hereafter 
SV98). These studies have produced some observational evi- 
dence in support of the variation of the IMF slope with metal- 
licity. Schaerer et al. (1999) reported the detection of WC 
stars in five WR galaxies, of which three have metallicities 
of Z = 0.004 and the rest two have Z = 0.008 and Z > 0.020 
respectively. By comparing the observations with the model 
predictions of SV98, these authors found that only the two 
objects with higher metallicities are consistent with the stan- 
dard model in which a Salpeter IMF is assumed, whereas the 
observed WR emission lines of the other three low metallic- 
ity objects can not be reproduced by this model and a flatter 
IMF is required in order to properly bring the models into 
agreement with the data. A study of three WR galaxies with 
Z = 0.004 by Huang et al. (1999) reached a similar conclu- 
sion that a top-heavy IMF with a = 1.00 is likely required 
to favor the formation of WR stars in these low-metallicity 
galaxies. Recently, Schaerer et al. (2000) investigated six 
metal-rich WR galaxies (Z ^ 0.020) and found that a very 
steep IMF with a ^ 3.30 is very unlikely for these objects. 
A more recent analysis of a larger sample of 14 WR galaxies 
by Fernandes et al. (2004) has concluded that there does exist 
some dependence of the IMF slope on galaxy metallicity. In 
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this study, low-metallicity galaxies are predicted to undergo a 
burst of star formation and show a fiat IMF varying between 
the top-heavy and the Salpeter form, whereas high metallicity 
galaxies either show a steeper IMF or experience an extended 
burst. 

The studies described above suggest that, at least for WR 
galaxies, the slope of the IMF is somewhat dependent on 
metallicity. However, the galaxy samples used in these pre- 
vious studies are not sufficiently large to give reliable con- 
straints on the relation between the IMF and the metallicity of 
galaxies. In this paper, we address this problem using a large 
sample of WR galaxies drawn from the Data Release Three 
(DR3; Abazajian et al. 2005) of the Sloan Digital Sky Survey 
(SDSS; York et al. 2000). SDSS is to date the most ambitious 
imaging and spectroscopic survey, providing a library of hun- 
dreds of thousands of galaxy spectra with a resolution of < 
3 A at 5000A and a median signal-to-noise ratio (S/N) of ~ 14 
per pixel. It thus has allowed one to select an unprecedented 
large, high quality spectral sample of WR galaxies, and sub- 
sequently to give reliable constrains on the relation between 
galaxy metallicity and the IMF slope. 

From DR3, we select a large sample of 174 galaxies with 
obvious WR features using a two-step method, in which we 
first select star-forming galaxies with evident He emission and 
then select WR galaxies from these candidates by visually ex- 
amining their spectra (§2). Next, for each WR galaxy we care- 
fully calculate its oxygen abundance and measure its WR and 
H(3 emission, as quantified by the equivalent width (EW) of 
the WR bump around 4650A, the EW of H/3 emission line, 
and the WRbump/H/3 intensity ratio (§3). We then investigate 
the metallicity dependence of the IMF slope by comparing 
these observed quantities with those predicted by evolution- 
ary synthesis models (§4). We discuss in §5 the possible bi- 
ases in our results caused by a few sources of systematic error. 
Finally, we summarize our results in the last section. 

2. SAMPLE SELECTION 
2.1. The SDSS spectroscopic sample 

The data analyzed in this study are drawn from the SDSS. 
The survey goals are to obtain photometry of one-quarter of 
the sky and spectra of nearly one million objects. Imaging 
is obtained in the u, g, r, i, z bands (Fukugita et al. 1996; 
Smith et al. 2002) with a special purpose drift scan camera 
(Gunn et al. 1998) mounted on the SDSS 2.5-meter telescope 
at Apache Point Observatory. The imaging data are pho- 
tometrically (Hogg et al. 2001) and astrometrically (Pier et 
al. 2003) calibrated, and used to select stars, galaxies, and 
quasars for follow-up fibre spectroscopy. Spectroscopic fi- 
bres are assigned to objects on the sky using an efficient tiling 
algorithm designed to optimize completeness (Blanton et al. 
2003). The details of the survey strategy can be found in York 
et al. (2000), and an overview of the data pipelines and prod- 
ucts is provided in the Early Data Release paper (Stoughton 
et al. 2002). More details on the photometric pipeline can be 
found in Lupton et al. (2001). 

Our parent sample for this study is composed of 374,767 
objects which have been spectroscopically confirmed as 
galaxies and have data publicly available in the SDSS DR3. 
It is a catalogue of local galaxies (mostly below z ~ 0.3 and 
a median z of 0.1), covering 4188 deg 2 on the sky. The spec- 
tra are obtained with two 320-fibre spectrographs that are also 
mounted on the SDSS 2.5-meter telescope. Fibers with 3 arc- 
sec in diameter are manually plugged into custom-drilled alu- 



minum plates mounted at the focal plane of the telescope. The 
spectra are exposed for 45 minutes or until a fiducial S/N is 
reached. The median S/N per pixel for galaxies in the main 
sample is ~ 14. The spectra cover an optical wavelength 
range from 3800 to 9200 A with an instrumental resolution 
of R = 1850-2200 (FWHM - 2.4 A at 5000 A ). 

2.2. Description of the selection method 

We aim to select a sample of galaxies showing evident WR 
features in their spectra (blue bump around 4650A and/or red 
bump around 5808A). An intuitively accessible way is to sys- 
tematically measure the EW of WR bumps and take those 
with large EWs as WR galaxies. However, this method will be 
extremely unefficient when dealing with large data sets, given 
that the WR galaxy population is just a very small fraction 
of the whole sample and in most cases systematical measure- 
ments of WR bump strength suffer from large uncertainties in 
the continuum determinations and/or from the contamination 
by the emission of non-WR stars. 

To avoid this difficulty, we instead opt for a two-step 
method in which we first select star-forming galaxies with 
evident He emission line as candidates, and then select WR 
galaxies by visually examining these candidates. There is no 
denying that this method will inevitably drop a number of 
WR galaxies, especially those with low S/N. However, this 
method can still work well for our purpose, because of the 
following two facts. First, as described in the introduction, 
WR galaxies are believed to be undergoing present or very 
recent star formation that produces massive stars. It is then 
natural to speculate that a sample of WR galaxies must be 
a subsample of star-forming galaxies and most WR galaxies 
must show high-order Balmer emission lines in their spectra 
because these lines are usually associated with massive stars 
that are expected to be numerous in WR galaxies. We use the 
He emission line rather than other Balmer lines such as HS 
in order to ensure that the high-excitation H II regions are se- 
lected. It should be pointed out that the He emission line is 
in fact a blend of three lines, including He, [Ne III] and He I 
lines. The latter two lines are often strong in active galac- 
tic nuclei (AGN). By using this blend we can find AGN with 
WR features in a harmonious way as normal galaxies, which 
are also interesting and will be discussed in a separated paper, 
at the cost of somewhat reducing our selection efficiency for 
the purpose of this paper. The second fact is that the parent 
sample is substantially large, allowing one to select a large 
enough sample of WR galaxies and thus to study their statisti- 
cal properties still with unprecedented accuracy, even though 
some objects must have been thrown away. 

2.3. Starlight removal and emission line measurements 

To select star-forming galaxies, we have performed a care- 
ful subtraction of the stellar absorption-line spectrum before 
measuring the nebular emission lines. Our methodology for 
modelling and subtracting the underlying starlight has been 
described in detail in Li et al. (2005) and in Lu et al. (2006). 
The method for measuring the emission-line parameters has 
been described in Dong et al. (2005). A set of six absorption- 
line templates with zero velocity dispersion are constructed in 
the first place, by applying the technique of Independent Com- 
ponent Analysis to the synthesis galaxy spectra from a newly 
released high-resolution evolutionary model by Bruzual & 
Chariot (2003). The physical meaning of the six templates 
could be easily understood through examining their spectra 
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visually (see Fig. 4 of Lu et al. 2006). The first template repre- 
sents the blue continuum of O stars as well as the Ca II H and 
Ca II K lines at 3933A and 3968 A. The second template is 
similar to the spectrum of B stars, while the absorption lines 
of neutral metals and molecules such as TiO are also identi- 
fied. The third template shows extremely strong Balmer ab- 
sorption lines, a Balmer jump at the blue end, and the Ca II 
triplet and more lines of neutral metals and molecules at the 
red end. The fourth and fifth templates are somewhat like hy- 
brids of F-K stars with stronger neutral metal and molecule 
lines. The last template is similar to the spectrum of M stars 
in the long-wavelength range but shows high-order Balmer 
absorption lines at short wavelengths. The spectral properties 
of the templates imply a tight correlation between the stellar 
population age and the templates (see Lu et al. 2006 for a de- 
tailed discussion). The spectra of all galaxies in the parent 
sample described above are then fitted with the six templates 
with iterative rejection of emission lines and bad pixels. Dur- 
ing the modeling, the template spectra are broadened to ve- 
locity dispersion from to 600 kms -1 using Gaussian kernel, 
and an average intrinsic reddening is added to the model by 
searching a range of color excess E(B-V) and assuming an 
extinction curve of Calzetti et al. (2000). For each spectrum, 
the best-fitting model is subtracted from the original spec- 
trum, yielding a pure emission-line spectrum from which the 
emission-line parameters (the central wavelength, line inten- 
sity and line width) are measured. 

As discussed in Li et al. (2005), our method of modeling 
and subtracting the underlying stellar component have several 
important advantages, and so the nebular emission lines can 
be accurately measured from the starlight-subtracted spectra. 
First, the templates are derived from a large sample with full 
coverage of spectral type. This ensures a very close match 
to the true underlying stellar population of different types of 
galaxies, including emission-line galaxies that contain much 
young stellar components. Second, the templates are not ob- 
tained directly from the observed spectra, but from the mod- 
eled one with zero velocity dispersion, and thus could well 
match the absorption-line profiles. Finally, all the emission 
lines with flux above 3er are carefully masked during the mod- 
eling. This is vital to the nebular line measurements, espe- 
cially for star forming galaxies. For each galaxy, the modeled 
spectrum is subtracted from the observed one, and the emis- 
sion lines are fitted with Gaussian functions (see Dong et al. 
2005 for details). Because emission lines and absorption lines 
are coupled, this procedure is performed iteratively. 

2.4. Selecting WR galaxies 

Based on the emission-line measurements, star-forming 
galaxies are selected from the subset of galaxies with the in- 
tegrated flux intensity F > 2 x 10~ 17 ergs~ l cm~ 2 in the four 
emission lines [OIII]A5007, H/3, [NIIJA6583 and Ha. Fol- 
lowing Kauffmann et al. (2003), a galaxy is defined to be a 
star-forming galaxy if 

log([OIII]/H/3) <0.61/(log([NII]/Ha)-0.05) + 1.3. (1) 
This gives rise to a total of 1 17,275 star-forming galaxies. 

The candidates of WR galaxies are then selected from 
these star-forming galaxies by requiring EW(He)> 5a, where 
EW(He) is the EW of He emission line and a is its error. As 
a result, a total of 7628 galaxies are selected as our candi- 
dates. As described above, the spectra have been well fitted 
and the best- fitting stellar component (continuum plus absorp- 
tion lines) has been separated from the emission-line compo- 
nent. Based on these fits, EW(He) and a can be measured 



in an easy and accurate manner. When measuring, we have 
adopted a rest-frame wavelength range of 3959A to 3979 A for 
He. For each object, we use the best-fitting stellar spectrum 
(see §2.2), but not the observed spectrum, to determine the 
continuum that will be used for calculating the EW, in order 
to avoid possible contamination by nearby nebular emission 
lines. We have adopted 3905-3915A for the nearby blue band 
and 4010-4030A for the red. In the computation, the flux fluc- 
tuations provided by the SDSS pipeline have been taken into 
account. 

Finally, we visually examine the spectra of the 7628 candi- 
dates and keep those with obvious WR features. This gives 
rise to a sample of 174 WR galaxies. In some cases, the se- 
lected WR galaxies show both the blue and red bumps, while 
in other cases only the blue bump is visible. We tabulate our 
sample in Table 1 so that the reader would be able to recover 
our results or perform other analyses using this sample. For 
each WR galaxy, we list in the table the spectroscopic ID 
(Column 1), the intensity ratio of various emission lines rel- 
ative to H/3 emission line (Columns 2-18), the EWs of H/3 
emission line (Column 20), the WRbump/H/3 intensity ratio 
(Column 21), the EWs of blue WR bump (Column 22) and 
the oxygen abundance (Column 23). Before these quantities 
being measured, the observed spectra have been corrected for 
the Galactic and the internal reddening (see the next section). 
The observed flux of H/3 emission line is listed in Column 19. 

3. OBSERVATIONAL MEASUREMENTS 

3.1. Correction of internal extinction 

In the previous section, we have described our method for 
subtracting the underlying starlight and for measuring the 
emission lines. Before modelling the stellar spectrum of 
galaxies, the observed spectra have already been corrected for 
the Galactic reddening, but not for the internal reddening. The 
internal reddening can be estimated using the Balmer decre- 
ments of the emission lines from H II regions. We didn't ap- 
ply this correction before the modelling, because our fitting 
method is aimed to model the stellar component of the galax- 
ies, but not the H II regions where the gas dominates the ob- 
served spectrum. It has been known that the mean color ex- 
cess estimated from the observed Balmer decrements is sys- 
tematically larger than the average of the stellar extinction 
(Chariot & Fall 2000; Li et al. 2005). Therefore, we choose 
to take into account the stellar reddening by a single average 
color excess during the modelling. In this way our program at- 
tempts to fit the observed uncorrected spectra (continuum plus 
absorption lines). After the modelling, we measure emission- 
line parameters, correct for intrinsic extinction using Balmer 
decrements, and use the corrected parameters for further anal- 
ysis. Our procedure of sample selection described above is 
not sensitive to this modification. However, when estimating 
metal abundance based on measurements of nebular emission 
lines, it is necessary to carefully correct for the internal red- 
dening. 

In practice, the internal reddening is usually corrected 
as a whole using the effective extinction curve t\ — 
t v (A/5500)" 07 , which was introduced by Chariot & Fall 
(2000). The effective V-band optical depth -ry, arising from 
attenuation by dust in a galaxy, can be written 

ln[F(Ha)/F(H^)] -ln[/(Hq)//(H/3)] 

TV [(A Hct /5500)- -MAtW5500)- - 7 ] ' 
where F(Ha)/F(H/3) is the ratio of the observed Ha emis- 
sion intensity relative to H/3, and /(Ha) //(H/3) is the theoret- 
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TABLE 1 
Parameters of the WR sample. 11 



MJD-Plale-Fibc(iOIIH10H9 H8 H7 H5Hi(OIIl]p III[P III]Ha[S HQS 11||N HON 110O IIDO II] 11/3 EWWREW 12+ method 
37237988368909791093401363 4959 5007656»71'B73 15548558373207330486 l<H0E>uim>65Ojg(O/H) 

(1) (2) a (3) (4) (5) (6) (7) (8) (9) (10) (1 1) (12)(13)(14)(15)(16) (17) (18)(19) b (20)(21)(22) (23) (24) 

52251-0751-629. ... 6.119.629.927.830.4 7.6 164.0515.3280.89.1 6.2 1.4 4.1 4.2 1.9 847.366.7 5.1 3.0 8.01 Te 

51900-0390-44591.42.6 5.3 17.828.724.817.7 3.6 151.2476.2285.47.913.6 7.2 20.7 2.3 1.91941.187.8 6.8 5.1 8.27 Te 

51900-0390-10606.9 101.39.221.013.2 ... 12.7 39.5286.O57.930.335.7102.0 1.1 0.5 274.318.422.13.5 9.04 N2 

52233-0753-094... 4.8 7.7 19.222.925. W6.8 2.4 132.9418.2287.0.4.710.6 4.4 12.7 2.0 1.7 2323.551.717.3 7.9 8.33 Te 

51817-0418-302... 4.16.4 1 8.226.225.016.8 5.6 154.4485.8282.515.611.4 3.0 8.7 2.0 1.5 2271.079.4 7.1 5.1 8.05 Te 

51913-0394-47293.33.9 6.116.728.724.945.311.1 172.9545.1278.48.2 6.1 1.3 3.7 1.1 1.0 723.886.8 2.7 2.1 7.80 Te 

51789-0398-294 6.420.022.225.147.3 1.1 94.8 297.6289.624.116.9 6.7 19.2 2.2 2.31462.476.116.311.0 8.42 Te 

51817-0399-32246.60.4 3.6 17.912.924.W6.2 0.5 29.5 92.0288. D7.829.730.085. 8 1.8 1.6 975.633.8 7.0 2.0 8.09 Te 

51882-0426-314 7.225.027.527. W8.5 10.0 120.3379.5 ... 9.3 7.2 0.9 2.6 1.4 0.7 989.8193.11.4 2.4 7.51 Te 

51900-0427-29225.20.12.4 16.912.422.543.3 1.6 64.4 201.6285.224.520.320.157.5 209.530.824.4 6.6 8.06 Te 



a) Columns 2-18 and 21 are intensity ratios relative to H/3 emission line, with /(H/3) — 100. Before these quantities being measured, the observed spectra have been corrected for the 
Galactic reddening, the internal reddening and the stellar absorption. See the text for details. 

b) The observed flux of H/3 emission line in unit of 10~ 17 erg s^cnT 2 , which has been corrected for the Galactic reddening, but not for the internal reddening. 

c) Dots denote the cases in which Ha or [O III] A A4959, 5007 lines are clipped, or the corresponding lines are not detected. 

d) Table 1 is presented in its entirety in the electronic edition of the Journal. A portion is shown here for guidance regarding its form and content. 



ical value of this ratio. It has been known that I(Ha)/I(H(3) 
does not keep constant, but varies with electron temperature 
T e , ranging from -2.7 at T e = 20000K to - 3.0 at T e = 
10000K. In order to take this variation into account, we de- 
termine the electron temperature T e using the observed inten- 
sity ratio [OIII]A4363/[0 III] AA4959, 5007, from which we 
then calculate I(Ha)/I(Hf3) using the method introduced by 
Storey & Hummer (1995). As a result, theoretical Ha-to-H/3 
ratios are successfully obtained for 152 WR galaxies whose 
T e can be determined. In general, the observed Ha/H/3 ra- 
tio is larger than the theoretical value, because H/3 emission 
is more heavily absorbed than Ha. However, out of the 152 
WR galaxies, there are 12 objects showing relatively smaller 
HafH/3 ratios than theoretical values. This could be due to 
clipped Ha lines or observational flux fluctuations. For these 
galaxies, we instead follow Izotov et al. (1994) to use H7/H/3 
intensity ratios to derive the extinction coefficient C(H(3) and 
the theoretical values of H7/H/3. In this way, 8 objects are 
corrected for the internal reddening. For the remaining 4 ob- 
jects whose observed H7/H/3 ratios larger than the theoretical 
values, we do not perform the correction. 

For the 22 galaxies without T e determinations, the correc- 
tion of the internal reddening is performed in the same way, 
except that the theoretical intensity ratios are set to typical val- 
ues of /(Ha/H/3) = 2.86 (18 objects) or /(H 7 /H/3) = 0.468 
(4 objects). 

Fig- HI compares the corrected hydrogen emission line ra- 
tios with the theoretical values for the galaxies that have T e 
determinations and have been performed the correction based 
on the Ha/H/3 intensity ratios. The crosses and filled circles 
show the observed and corrected ratios respectively. The lines 
are for the theoretical values. As can be seen, the intensity 
ratios of hydrogen lines Ha, H7 and HS relative to H/3 are 
corrected substantially and change to be consistent with the 
theoretical values. 

Before going to the next section, it is necessary to 
find out the objects in our sample with clipped Ha or 
[O III]AA4959,5007 lines in their spectra, because the oxy- 
gen abundance of galaxies will be determined based on 
these lines. By comparing the observed and/or corrected 
emission line ratios with the corresponding theoretical val- 
ues, we find that 2 objects have been clipped for Ha 
line, 1 for [OIII]AA4959,5007, and 2 for both Ha and 
[OIII]AA4959,5007 lines. 
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FIG. 1 . — Intensity ratio of hydrogen emission lines as a function of elec- 
tron temperature T e , for WR galaxies in our sample that have been corrected 
for internal reddening. The observed ratios are plotted as crosses. The filled 
circles show the results that have been corrected for internal reddening using 
the method described in the text. The curves are for the theoretical values. 
[See the electronic edition of the Journal for the color version of this figure.] 



3.2. Oxygen Abundance 

There are several empirical methods that have been pro- 
posed to estimate the oxygen abundance of galaxies. In this 
paper we adopt two of them: the classic T e method and the 
N2 method. In another paper (Shi et al. 2005), we have ap- 
plied these methods to a sample of blue compact galaxies and 
discussed in detail the possible systematic differences among 
different methods. We present below a brief introduction to 
these methods and refer the reader to Shi et al. (2005) for a 
detailed description. 

In the T e method, one first determines the elec- 
tron temperature T e and density n e of the gas us- 
ing intensity ratios [OIII]A4363/[0 III] AA4959, 5007 and 
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[S II]A6716/[S II] A673 1 respectively, and then computes the 
abundances of + and O" 1 " 2 using the IRAF package NEBULAR 
(Shaw & Dufour 1995). In the computation, the + abun- 
dance is determined by the intensity of [OII]A3727 line or 
[OIII]AA7320,7331, while the abundance of +2 is deter- 
mined by two [O III] lines centered at 4959 and 5007A. More- 
over, the fact that He IIA4686 line is usually detected in the 
spectra of WR galaxies indicates that part of the oxygen in 
these systems must be in O" 1 " 3 stage. Following Skillman & 
Kennicutt (1993), we estimate the O" 1 " 3 abundance by assum- 
ing +3 /(0 + + +2 ) = He +2 /He + , where the helium abun- 
dances are determined using the He IA6678 and He IIA4686 
lines (Clegg 1987; Skillman & Kennicutt 1993). Therefore, 
the total abundance of oxygen in WR galaxies is derived by 
O/H = (0 + + +2 + +3 )/H + (Izotov et al. 1994; Skillman & 
Kennicutt 1993). 

As mentioned above, we have determined the electronic 
temperature T e for 152 out of 174 galaxies in our sample. We 
thus determine the oxygen abundance using the T e method for 
these galaxies. In this procedure, the O" 1 " abundance is deter- 
mined by the intensity of [OII]A3727 for 77 objects, while 
for the other 75 with the non-observed [OII]A3727 line, the 
[OIII]AA7320,7331 lines are used. 

For the remaining 22 galaxies without T e determinations, 
we have to turn to a second method, because some of the 
emission lines mentioned above, which are necessary in the 
T e method, cannot be measured from their spectra due to the 
limited wavelength coverage (e.g. in case of [OII]A3727), 
the low S/N (e.g. in case of [OIII]A4363) or clipping of 
[OIII]AA4959,5007. For these galaxies, we adopt the so- 
called N2 method (van Zee et al. 1998), in which the oxygen 
abundance is estimated by 

12 + log(0/H) = 1.021og[/([NII])//(Ha)]+9.36. (3) 

Where /([Nil]) is /([NII])(A6548 + A6583). It should be 
pointed out that for 3 objects whose Ha lines are clipped, 
the Ha fluxes are fixed to be 2.86x/(H/3). We list in Table 
1 the oxygen abundance for each object (Column 23) and the 
method adopted for it (Column 24). 

As pointed out by Shi et al. (2005), the oxygen abundances 
derived from the N2 method are systematically larger by ~ 
0.2dex than those from the T e method. We will show that this 
systematic difference does not affect our conclusion on the 
metallicity dependence of the IMF slope late. 

3.3. The WR and Hf3 emission 

In this paper we make comparisons with models for three 
quantities: the EW of H/3 emission, the EW of WR bump 
around 4650A (the blue bump), and the blue bump-to-H/3 in- 
tensity ratio WRbump/H/3. As first pointed out by Copetti, 
Pastoriz, & Dottori (1986), the EWs of hydrogen recombina- 
tion lines, particularly H/3, can be used as an accurate indica- 
tor of the age of the population in the case of an instantaneous 
burst of star formation. In addition, as pointed out by SV98, 
EW(H/3) as a function of metallicity can be combined with 
WR lines/H/3 intensity ratios to determine the WR and O star 
populations. For the WR features, we concentrate on the blue 
bump only because this feature is visible in all the spectra (see 
§2.2). 

We have described in detail our method of measuring H/3 
flux in §2.2 and the method of correcting for the effect of 
internal extinction in §3.1. The EW of H/3 is measured 
in a similar way. However, The situation becomes difficult 



when one trying to systematically measure the EW and in- 
tensity of WR bumps. As already mentioned in §2.2, sys- 
tematical measurements of WR bump strength usually suf- 
fer from large uncertainties in the continuum determinations 
and/or from the contamination by the emission of non-WR 
stars such as [FeIII]A4658, HeIIA4686, He I+[ArIV]A471 1, 
and [ArIV]A4740. We thus choose to manually measure the 
WR bump for all our WR galaxies. First, for each object, 
we use the best-fitting stellar spectrum, but not the observed 
spectrum, to determine the continuum. The continuum is then 
determined by fitting a cubic spline to seven wavelength chan- 
nels (4020, 4510, 5313, 5500, 6080, 6630, 7043A), which are 
thought to be fully free of nebular and stellar lines (Kong et 
al. 2002). Next, we fit the non-WR emission lines mentioned 
above with Gaussian functions and subtract them from the ob- 
served spectrum. We then measure the EW and the integrated 
flux of the blue bump with rest-frame wavelength limited to 
4600-4750A. 

4. COMPARISON WITH MODEL PREDICTIONS 

In this section, we compare the observed WR features and 
H/3 emission with predictions by evolutionary synthesis mod- 
els of SV98. In the comparison, we divide our sample into 
four subsamples according to oxygen abundance and inves- 
tigate whether/how the slope of the IMF depends on galaxy 
metallicity. We first briefly describe the models, and then 
present the results of the comparison. 

4. 1 . Predictions by Evolutionary Synthesis Models 

The model predictions used in this paper are provided by 
S V98. These authors constructed evolutionary synthesis mod- 
els for young starbursts, using stellar evolution models, theo- 
retical stellar spectra, and a complication of observed emis- 
sion line strengths from WR stars. They explicitly distinguish 
between the various WR subtypes whose relative frequency 
is a strong function of metallicity. The observational features 
predicted by their models allow a detailed quantitative deter- 
mination of the WR star population in a starburst region from 
its integrated spectrum and provide a means of deriving the 
burst properties (e.g. duration and age) and the slope of the 
IMF of young starbursts. The model predictions should pro- 
vide the most reliable determinations to date, and can be used 
to investigate the variation in stellar properties with metallic- 
ity. 

In the models, the input parameters include the burst du- 
ration Af, the metallicity Z, the IMF, and the total mass of 
stars formed in the burst. As pointed out by SV98, the last 
parameter serves only as a normalization constant and is not 
relevant for the quantities being studied. For the IMF, a power 
law dN I dM °c M~ a between the upper and low cutoff masses, 
M up and M/ ow respectively, is adopted. Following SV98, we 
set the cutoff masses to M up = 12OM and M low = O.8M . 
For the slope of the IMF, we consider three different cases: a 
Salpeter slope with a = 2.35, a steeper slope with a = 3.30, 
and a flatter slope with a = 1 .00. In each case, we assume an 
instantaneous burst with Af = 0. In order to test the possible 
dependence on burst duration, we also consider in the case 
of Salpeter IMF four extended burst models with durations 
Af =1, 2, 4 and 10 Myr. Finally, in order to study the depen- 
dence on metallicity, we consider four cases of metallicity: 
Z=0.001, 0.004, 0.008, and 0.02. 

With each set of input parameters, the models predict many 
important quantities which can be easily compared with ob- 
servations. These quantities are classified into three groups: 
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population statistics (the number of stars with different spec- 
tral types), nebular quantities (the strength of nebular recom- 
bination lines: Ha, H/3, HeIIA4471, and HeIIA4686), and 
WR emission lines. In this paper we study three of them: 
EW(H/3) (the EW of H/3 emission line), EW(4650) (the EW 
of the blue WR bump), and WRbump/H/3 (the ratio of blue 
WR bump intensity relative to H/3). 

4.2. Metallicity dependence of the IMF 

In order to study the metallicity dependence of the IMF and 
to make comparisons with the model predictions described 
above, we divide our sample into four subsamples according 
to oxygen abundance. Following Guseva et al. (2000), we as- 
sume the solar oxygen abundance [12 + log(O/H)] = 8.93 
and use the oxygen abundance ranges listed below for the 
sample division. 



1. 


Z= 


=0.001: 


12 + log(0/H) <7.93; 




2. 


Z= 


=0.004: 


7.93<12 + log(0/H)<! 


3.43; 


3. 


Z= 


=0.008: 


8.43<12 + log(0/H)<! 


3.63; 


4. 


z= 


=0.020: 


12 + log(0/H) > 8.63. 





Figure |2] shows the average spectrum of WR galaxies in 
each of the four metallicity subsamples. In the calculation, 
each observed spectrum is weighted by its S/N and normal- 
ized by the flux at 5000A. Panels from top to bottom corre- 
spond to different metallicities, as indicated in the left-hand 
panels. The number of objects included in each subsamples is 
also indicated. The panels at left-hand show the results ob- 
tained using the observed spectra that have been corrected 
only for Galactic reddening, while the right-hand panels are 
for those also corrected for the internal reddening. The inset in 
each panel shows the blue WR bump around 4650A. It can be 
seen that, if the internal extinction was not taken into account, 
low metallicity population shows harder continuum, as well 
as stronger nebular emission lines, than does high metallicity 
population. After the internal reddening being corrected for, 
the difference in the shape of the continuum disappears and 
the continuum slope changes to be very close to each other. 
This is reasonable because the dust content in high metallic- 
ity galaxies is expected to be larger than that in low metallicity 
ones. As a result, dust extinction will modify the continuum 
more significantly for high metallicities. It is also interesting 
to see that, the decrease of the strength of nebular emission 
lines does exist for large metallicities, regardless of whether 
the internal reddening is taken into account or not. 

These observational results are shown more clearly in Fig- 
ure|3] where we have plotted three quantities (from top to bot- 
tom: EW(4650), EW(H/3) and WRbump/H/3) as a function of 
oxygen abundance for all our WR galaxies (crosses). As ex- 
pected, the observed WRbump/H/3 intensity ratio decreases 
with decreasing metallicity. This result was first found by 
Arnault, Kunth, & Schild (1989) in an attempt to quantify 
the constrains on the WR populations in starbursts, and were 
largely confirmed by later studies. The result is believed to be 
due to the strong influence of metallicity on stellar mass loss, 
which leads to a significant decrease in the WR population 
at low metallicities (Maeder, Lequeux, & Azzopardi 1980). 
Moreover, the rapid decrease of EW(H/3) with increasing oxy- 
gen abundance is clearly seen from this figure. In contrast, we 
see no significant correlation between EW(4650) and metal- 
licity. The increase of WRbump/H/3 with metallicity is thus 



in large part due to the rapid decrease of H/3 for large metal- 
licities, which is well consistent with theoretical expectations 
(see below; also see Arnault et al. 1989). 

Figure [5] also compares these quantities to predictions 
by the three instaneous bust models with different IMFs 
(a =1.00, 2.35, 3.30) and metallicities (Z =0.001,0.004, 
0.008, 0.020). We find that models with flatter IMFs tend to 
predict larger values for all these quantities and at all metal- 
licities. This is reasonable because a flat IMF implies a rel- 
atively large population of massive stars, leading to strong 
WR and nebular emission. In particular, such metallicity de- 
pendence is much stronger in EW(4650) and WRbump/H/3 
than in EW(H/3). We also see that the observed trends in 
the three quantities with oxygen abundance can generally be 
reproduced by all the models. However, the strength of the 
predicted trends is different for models with different IMFs. 
For example, a flatter IMF will lead to more rapid increase 
of WRbump/H/3 with metallicity. As a result, the difference 
between model predictions becomes much more remarkable 
at high metallicities, reaching a factor of ~ 5 in EW(4650) 
and ~ 3 in WRbump/H/3 on average at Z > Z . It is thus 
inevitably to see that a model with a fixed IMF can not simul- 
taneously match the observations at all metallicities. At the 
highest metallicity (Z = 0.02), only the model with the steep- 
est IMF (a = 3.30) can well match the observed EW(4650) 
and WRbump/H/3. In case of Z = 0.001, however, just the 
reverse is likely to be true. In another word, Fig.[3]indicates 
that a dependence of the IMF slope on galaxy metallicity is 
required to explain the observed trends in WR and H/3 emis- 
sion. 

Figures 0] and |5] (black points in the top panels) show 
the distribution in spaces of EW(4650) versus EW(H/3) and 
WRbump/H/3 versus EW(H/3) respectively, for WR galaxies 
in different metallicity subsamples (as indicated above each 
panel). In each panel, the average location of the objects is 
plotted as a green point and the 1 a dispersion between them 
is indicated by the error bars. The solid lines show the predic- 
tions by the three instaneous burst models, with red, magenta 
and blue colours respectively corresponding to models with 
IMF slope a=1.00, 2.35 and 3.30. The magenta lines in other 
styles are for the four extended burst models with fixed IMF 
slope (a = 2.35), but various burst durations: Af= 1 (dotted), 
2 (dashed), 4 (long dashed), 10 Myr (dotted-dashed). A clear 
dependence of the IMF slope on metallicity can be seen from 
the two figures, in that WR galaxies with higher metallici- 
ties exhibit steeper IMF. As can be seen, most objects with 
Z = 0.02 are located well below the a = 3.30 curve, implying 
that these galaxies have IMFs similar to or even steeper than 
the steepest IMF that we have probed. In contrast, the objects 
of Z = 0.001 are dominantly above the a = 1.00 curve, and 
so should be associated with very flat IMFs. At intermedi- 
ate metallicities, the scatter in the observational distributions 
is very large, but the observations are still on average consis- 
tent with the Salpeter IMFs, the ones having the intermediate 
slope. 

5. SOURCES OF SYSTEMATIC ERROR 

In the previous section, we have seen a clear dependence of 
the IMF slope on metallicity by comparing our measurements 
of WR features and H/3 with the model predictions provided 
by SV98. However, there are possible uncertainties in this 
result, caused by a few sources of systematic error either in 
the models or in the observational measurements. We discuss 
these in detail in this section. 
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FIG. 2. — Average spectrum of WR galaxies in different metallicity subsamples. Panels from top to bottom correspond to different metallicities, as indicated in 
the left-hand panels. The number of objects in each subsample is also indicated. The panels at left-hand are for the results obtained from the observed spectra that 
have been corrected only for Galactic reddening, whereas the right-hand panels are for those also corrected for the internal reddening. The inset in each panels 
shows the WR bump around 4650A. 



5.1. Aperture effects 

Each of the SDSS spectra was obtained using a 3 arcsecond 
diameter fibre that had been positioned as close as possible 
to the centre of the object being studied. One may thus won- 
der to what extent our measurements of EW(4650), EW(H/3) 
and WRbump/H/3 may be biased because the measurements 
are not accurately reflecting the stellar and gaseous content of 
the whole starburst region. For example, because of aperture 
effects, both EW(4650) and EW(H/3) may be somewhat un- 
derestimated. As a result, the locations of galaxies on each 



panel of Figure|4]may be biased toward lower-left corner and 
thus mismatch the theoretical curves. 

In order to address the effect of aperture bias, we define a 
subsample of "safe" objects as those fully enclosed within the 
spectroscopic aperture. We use the Petrosian radii in the five 
SDSS photometric bands provided by the SDSS pipeline to 
quantify the angular size and select "safe" objects by requir- 
ing that at least one of the five radii is smaller than the radius 
of the aperture. This results in a subsample of 23 WR galax- 
ies whose spectra are not suffering from the aperture effect. 
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FIG. 3. — Equivalent width (EW) of WR bump around 4650A (top), EW of H/3 emission line (median), and the intensity ratio of WR bump relative to H/3 
(bottom) are plotted as a function of oxygen abundance. The crosses are for all the WR galaxies in our sample, while the filled circles are for "safe" objects, 
defined as those with at least one of the five Petrosian radii smaller than the radius of the SDSS spectroscopic aperture (1.5"). The symbols with error bars show 
the average values with \-u variance, which are predicted by models with different IMF slopes, as indicated in the top panel. Dashed vertical lines denote the 
oxygen abundance ranges and the model predictions are plotted at the median of each range. For clarity, triangles, squares and pentagons are slightly shifted 
along the horizontal axis. 



Though much smaller in size, this subsample will provide a 
critical test on the results obtained from the whole sample. 

In Figure|3l the "safe" objects are highlighted as filled cir- 
cles. In case of EW(4650), it is likely that there is no sig- 
nificant bias in the whole sample relative to the "safe" sam- 
ple. In contrast, a fraction of the whole sample shows much 
smaller EW(H/3) compared to the "safe" sample. We thus 
conclude that all emission of WR stars is likely inside the fi- 



bre, whereas aperture bias could well be a problem for H/3. 
This can be understood as H/3 emission is more extended. As 
a result, WRbump/H/3 intensity ratio is overestimated to vary- 
ing degrees for some of the objects. 

The relations between the three quantities obtained from the 
"safe" sample are plotted in the second rows of panels of Fig- 
ures[4]and [5] The symbols and the lines are the same as in top 
panels and the number of objects included in each metallicity 
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FIG. 4. — Distribution in the space of EW(4650) versus EW(H/3) for WR galaxies in different metallicity subsamples (black points). Panels from left to 
right correspond to different metallicities, as indicated above the top panels. The top panels are for the whole sample. The median panels are for objects with at 
least one of the five Petrosian radii smaller than the radius of the SDSS spectroscopic aperture. The bottom panels are for objects with their oxygen abundances 
determined from T e method. In each panel, the average location of the objects is plotted as a green point and the 1 — er dispersion between them is indicated by the 
error bars. The solid lines show the predictions by the three instaneous burst models, with red, magenta and blue colours, respectively, corresponding to models 
with IMF slope a=l.()0, 2.35 and 3.30. The magenta lines in other styles are for the four extended burst models with fixed IMF slope (a = 2.35) but various burst 
durations: At- 1 (dotted), 2 (dashed), 4 (long dashed), 10 Myr (dotted-dashed). 



subsample is also indicated. It is encouraging to see that the 
metallicity dependence of the IMF slope still exists. We are 
thus convinced that the effects of aperture do not change our 
main results. 

5.2. Oxygen abundance 

In our sample, the oxygen abundance of 152 objects is de- 
termined from the T e method, while for the others we have 
adopted the N2 method. As mentioned above, the measure- 
ments based on the N2 method are systematically larger than 
those based on the T e method, by a factor of ~ 0.2 dex. More- 
over, for 27 galaxies with [O III]A4959/H/3< 0.7, the EW of 
H/3 is generally low and the [O III]A4363 line is weak, hence 
their oxygen abundances might be underestimated (Izotov et 
al. 2004). It is thus necessary to test how reliable the methods 
are and whether and how they affect our main results. 

For T e method, it is possible to address this problem in a sta- 
tistical way by comparing the Ne/O, S/O and Ar/O abundance 
ratios and the average values found for galaxies with different 
oxygen abundances. If these ratios are not consistent with the 
average values, the metallicity might be somewhat underesti- 
mated (e.g. Izotov et al. 2004). Adopting a two-zone pho- 
toionized H II region model described by Izotov, Thuan, & 



Lipovetsky (1994), for each of the 152 objects whose oxygen 
abundances are determined from the T e method, we calculate 
its Ne ++ /0 ++ , S + and S" 1 " 1 " using expressions from Pagel et al. 
(1992), and Ar ++ , Ar +++ using IRAF NEBULAR package (see 
Shaw & Dufour 1995). We then calculate their Ne/O, S/O and 
Ar/O abundance ratios, using the ionization correction factors 
(ICFs) presented by Izotov, Thuan, & Lipovetsky (1994). The 
results are shown in Figure |6] The dotted lines are for the av- 
erage values found by Izotov & Thuan (1999). It is seen that 
the abundance ratios are well consistent with the average val- 
ues at all oxygen abundances. 

Figure[6]demonstrates that our calculation of oxygen abun- 
dance based on the T e method is reliable and robust. We 
thus select the 152 objects with oxygen abundance determined 
from the T e method to construct a second "safe" sample, in or- 
der to exclude any possible biases in our results caused by the 
N2 method. Similarly, we plot the results for this "safe" sam- 
ple in the third rows of panels of Figures|4]and |5]and obtain 
a conclusion that our results are not significantly biased by 
the uncertainties in oxygen abundance determinations. Note 
that there are no objects included in the Z = 0.020 subsam- 
ple. This is due to the fact that the application of T e method 
requires an accurate determination of the [O III]A4363 inten- 
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FIG. 5. — Distribution in the space of WRbump/H/3 versus EW(H/3), for WR galaxies in different metallicity subsamples. The symbols and the lines are the 
same as in Figurel4l 
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FIG. 6. — The Ne/O, S/O, and Ar/O abundance ratios of 152 WR galaxies 
whose oxygen abundance are determined from the T e method. 



sity, but in general this line is very weak or invisible in the 
spectra of high-metallicity galaxies. Nevertheless, the rest 
three panels still show marked increase of IMF slope with in- 
creasing metallicity, and thus our conclusion about the metal- 
licity dependence keeps unchanged. 

5.3. Spectral modelling 

In this paper, the EWs are measured adopting the best- 
fitting stellar continuum. However, the contribution of the 
ionized gas emission could be high in H II regions, especially 
for those with high EW(H/3). The most recent paper by Gu- 
seva et al. (2006), where the ionized gas emission is consid- 
ered in the fitting procedure, shows that the flux fraction of the 
gaseous continuum near H/3 can reach ~ 50% in H II regions 
with EW(H/3)~ 400A. Therefore, the ionized gas continuum 
template should be included in the spectral modelling, other- 
wise the EWs will be underestimated. Since our templates 
are constructed from simple stellar populations, our method 
currently can not be able to include the gas continuum in the 
fitting. Moreover, the absolute majority of our sample are low 
redshift objects in which the Balmer continuum is not within 
the SDSS wavelength coverage. It is thus very difficult to 
break the degeneracy between gaseous continuum and very 
young stellar populations only using the optical spectra with- 
out this feature. 

In addition, the measured EWs are compared with the 
model predictions for instantaneous or short bursts, and thus 
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the EWs should be obtained adopting the continuum level 
only of the youngest template that is most likely associated 
with the WR stars. Therefore, adopting the best-fitting stellar 
continuum will also lead to underestimated EWs. If the ratio 
of the youngest template continuum to the total one systemat- 
ically varies at different metallicities, then the results shown 
in Figures and [5] may also be biased. We have checked 
the relative contributions to the best-fitting spectra by the first 
two templates, which roughly represent the spectra of O and 
B stars respectively (see §2.3), and found no tendency of these 
contributions with metallicity. Therefore, although the EWs 
have somewhat been underestimated, the overall correlation 
between the measured EWs and the metallicities can still ex- 
ist and thus the trend of the IMF slope with metallicity can 
still be true. 

5.4. Synthesis models 

In this paper, the theoretical predictions are made using the 
evolutionary synthesis models of SV98. It is worthy noticing 
that their models are based on the evolution of non-rotating 
stars. Stellar rotation is expected to predict longer WR stage, 
lower WR star masses, and thus larger WR populations. As 
mentioned above, the strong influence of metallicity on stel- 
lar mass loss leads to a significant decrease in the WR pop- 
ulation at low metallicities, resulting in rapid decrease of 
WRbump/H/3 intensity ratio at these metallicities. Therefore, 
if rotating stars are included in the models, both EW(4650) 
and WRbump/H/3 would be larger than that predicted by the 
current models. As a result, the predicted curves in the low- 
metallicity panels of Figures 0] and [5] would go upward, and 
the observed decrease of the IMF slope at low metallicities 
would be less pronounced. New models of rotating stars are 
thus needed to be worked out, in order to provide more reli- 
able predictions, especially for low-metallicity populations. 

Moreover, S V98 have compiled average line fluxes for WR 
stars in the LMC (for models with Z < Z©) or the Galaxy (for 
models with Z > Z Q ), which represent the standard calibration 
for WR populations in external galaxies. There has, however, 
been a recent study of the luminosities of WR emission lines, 
reaching the indication that WR stars at SMC metallicities 
possess lower optical line luminosities than those in the LMC 
(Crowther & Hadfield 2006). This could be another source of 
systematic error which may also leads to underpredicted WR 
populations. 

It is thus important to answer the questions that to what 
extent the two effects described above would change the pre- 
dictions at low metallicities and which effect should be more 
dominant in such bias. However, at the moment, there are no 
published evolutionary synthesis models with rotating stars 
and more recent WR line luminosities. Nevertheless, one can 
still expect that these effects do not significantly change the 
results at high metallicities, for which the observed metallic- 
ity dependence of the IMF slope can still exist. 

6. SUMMARY 

In this paper, we have analyzed the WR features and H/3 
emission of 174 WR galaxies drawn from the Sloan Digital 
Sky Survey. By comparing these observed quantities with 
predictions by evolutionary synthesis models of Schaerer & 
Vacca (1998), we investigate whether and how the stellar ini- 
tial mass function depends on metallicity. Our sample is con- 
structed using a two-step method, in which we first select 
candidates as the star-forming galaxies with evident He emis- 
sion lines and then select WR galaxies by visually examine 



the spectra of the candidates. Although this method has in- 
evitably exclude a number of objects (mainly with low S/N), 
it is still work well for our purpose and gives rise to a large 
enough, high quality sample of WR galaxies. We have care- 
fully examine the possible sources of systematic error either 
in models or in our observational measurements and shown 
that these sources do not change our main results. 

The conclusions of this paper can be summarized as fol- 
lows. 

1 . The continuum of the average spectrum of WR galaxies 
shows no significant dependence on metallicity. In con- 
trast, low-metallicity galaxies exhibit in their spectra 
much stronger nebular emission than high-metallicity 
galaxies. 

2. The observed equivalent width (EW) of H/3 decreases 
with increasing metallicity, whereas there is no signifi- 
cant trend in EW of WR blue bump around 4650A with 
metallicity. 

3. We have found a clear dependence of the IMF slope on 
metallicity in that the IMF slope increases with increas- 
ing metallicity. Such dependence has been probed by 
many previous studies, but it is much more pronounced 
in this paper due to the improvement in the observa- 
tional data. 

It should be pointed out that there is still room for improve- 
ment of this study. First, the method of sample selection could 
be developed as to fully take advantage of the large number 
of SDSS galaxies. Second, more accurate measurements of 
WR lines and new models incorporating stellar rotation are 
needed to refine the results, especially for low metallicities. 
Third, the aperture effects should be corrected in a statisti- 
cally robust way incorporating the stellar contents obtained 
in our starlight fitting procedure and the method of spectral 
modelling could be improved. Finally, more quantities such 
as the relative numbers of WR stars with different subtypes 
could be included in the comparison to give more subtle and 
convincing results. These will be the aims of subsequent pa- 
pers. 
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